virology 177, 646-656 (1990) 


Establishing a Genetic Recombination Map for Murine Coronavirus 

Strain A59 Complementation Groups 

RALPH S. BARIC,*' 1 KAISONG FU,* MARY C. SCHAAD,* and STEPHEN A. STOHLMANt 


* Department of Parasitology and Lab Practice, School of Public Health, University of North Carolina at Chapel Hill, Chapel Hill, 
North Carolina 27599-7400; and \ Department of Neurology and Microbiology, University of Southern California, 

School of Medicine, 2025 Zonal Avenue, Los Angeles, California 90007 

Received December 7, 1989; accepted April 20, 1990 

MHV-A59 temperature-sensitive mutants, representing one RNA + and five RNA~ complementation groups, were 
isolated and characterized by genetic recombination techniques. Maximum recombination frequencies occurred under 
multiplicities of infection greater than 10 each in which 99.99% of the cells were co-infected. Recombination frequen¬ 
cies between different ts mutants increased steadily during infection and peaked late in the virus growth cycle. These 
data suggest that recombination is a late event in the virus replication cycle. Recombination frequencies were also 
found to range from 63 to 20,000 times higher than the sum of the spontaneous reversion frequencies of each ts 
mutant used in the cross. Utilizing standard genetic recombination techniques, the five RNA - complementation groups 
of MHV-A59 were arranged into an additive, linear, genetic map located at the 5' end of the genome in the 23-kb 
polymerase region. These data indicate that at least five distinct functions are encoded in the MHV polymerase region 
which function in virus transcription. Moreover, using well-characterized ts mutants the recombination frequency for 
the entire 32-kb MHV genome was found to approach 25% or more. This is the highest recombination frequency 
described for a nonsegmented, linear, plus-polarity RNA virus. © 1990 Academic Press, inc. 


INTRODUCTION 

Temperature-sensitive (ts) mutants of animal viruses 
have been used to elucidate many of the basic mecha¬ 
nisms of virus transcription, replication, assembly, and 
release in eukaryotic cells. Not surprisingly, ts mutants 
have been described for most oncogenic and nonon- 
cogenic RNA and DNA viruses (Fenner, 1970; Ghen- 
don, 1972). The gene order for several DNA-containing 
viruses has been established by DNA recombination 
analysis utilizing ts mutants and standard genetic tech¬ 
niques (Epstein era/., 1963; Ghendon, 1972; Ritchie, 
1973). Among those RNA viruses possessing seg¬ 
mented genomes, intercistronic reassortment be¬ 
tween wild-type or ts genome segments has identified 
the location and function of many reovirus, bunyavirus, 
and influenza virus genes (Fields, 1981). Reassortment 
frequencies for the influenza viruses and reoviruses 
vary from 2 to 42% depending on the virus strain and 
genome segment (Fields and Joklik, 1969; Fields, 
1981; Scholtissek etal., 1978). 

Evidence for homologous RNA recombination 
among RNA viruses with nonsegmented genomes has 
been reported for polioviruses, apthoviruses, and coro- 
naviruses (Cooper, 1968; Lai era/., 1985; Lake era/., 
1975; Mackenzie et at., 1975). Recombination fre¬ 
quencies between poliovirus ts mutants approach 
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2.2% or roughly 1% per 1700 nucleotide pairs of dou¬ 
ble-stranded RNA (Cooper 1968, 1977; Cooper era/., 
1975). Poliovirus recombination probably occurs dur¬ 
ing negative strand RNA synthesis by a copy choice 
mechanism (Kirkegaard and Baltimore, 1986); how¬ 
ever, the apthovirus and coronavirus recombination 
mechanisms are unclear. Numerous RNA recombinant 
viruses of mouse hepatitis virus (MHV), a member of 
the Coronaviridae, have been isolated and character¬ 
ized by T1 fingerprint analysis or sequencing of cross¬ 
over sites (Keck et a!., 1987, 1988b; Lai et at., 1985). 
Among different MHV strains, a region of polymor¬ 
phism and deletion has been detected in the MHV S 
glycoprotein gene (Parker et a!., 1989). A high fre¬ 
quency of RNA recombination has been suggested 
previously, but an exact frequency for the entire MHV 
genome has not been determined (Makino et at., 
1986b). Analysis of MHV-infected cells has revealed 
the presence of discrete RNA products bound to and 
dissociated from the replicative intermediate RNA 
which might represent transcription and RNA recombi¬ 
nation intermediates (Baric et a!., 1983, 1985, 1987). 
The precise role of these intermediates in transcription 
and recombination has yet to be determined. Se¬ 
quence analysis also suggests that nonhomologous 
recombination between wild-type MHV and influenza 
C virus may occur during mixed infection (Luytjes etal., 
1988). The mechanism for nonhomologous recombi¬ 
nation is unclear. 
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Recombination between the genomes of nonseg- 
mented RNA viruses provides a convenient tool to map 
the genetic loci of individual ts lesions and viral genes. 
Utilizing standard genetic recombination techniques, 
poliovirus and apthovirus ts mutants have been ar¬ 
ranged into an additive, linear, genetic map with mu¬ 
tants at different locations differing in physiologic func¬ 
tion (Cooper et at., 1975; Lake et ai, 1975). This has 
permitted the localization of many structural and non- 
structural gene functions to particular proteins. Tem¬ 
perature sensitive mutants of MHV have been isolated 
and characterized by complementation analysis and 
RNA phenotype (Koolen et ai, 1983; Leibowitz et at., 
1982; Martin ef a/., 1988; Schaad eta!., 1990). Utilizing 
the ts mutants described in the accompanying manu¬ 
script (Schaad et ai, 1990) and standard genetic re¬ 
combination techniques, we have formulated the first 
additive, linear, genetic map for MHV-A59 ts mutants. 
These data suggest that five or more genetic functions 
are encoded at the 5' end of the genome which partici¬ 
pate in regulating RNA synthesis. Moreover, utilizing 
well-characterized ts mutants from different comple¬ 
mentation groups, our data suggests that the actual 
recombination frequency for the MHV genome may ap¬ 
proach 25%. This is the highest homologous RNA re¬ 
combination frequency reported for a nonsegmented 
RNA virus. 


METHODS 

Virus propagation and assay procedures 

Temperature sensitive mutants representing one 
RNA + and five RNA” complementation groups of MHV- 
A59 were used throughout the course of this study. 
Virus stocks were propagated at 32° in 150-cm 2 flasks 
containing DBT cells as described in the accompany¬ 
ing paper (Schaad et at., 1990). Conditional lethal mu¬ 
tants used in this study represent the RNA - comple¬ 
mentation group A (LA3, LA6, LAI6, NC8); group B 
(LAI6, NC2, NC11); group C (LA8, LA9, NCI, NC10); 
group D (LA10); and group E (LAI8, NC4, NCI 2). The 
group F RNA + mutants used in this study include LA7 
and LAI 2. Plaque assays were performed at 32 or 
39.5° in Dulbecco’s modified essential medium 
(DMEM) (Sigma) in DBT cells containing 10% Nu-se- 
rum (Collaborative Research Inc.), 1% antibiotic/anti¬ 
mycotic (GIBCO), and 0.8% agarose. All plaque assays 
were quantified by staining with neutral red for 2 hr at 
28-36 hr postinfection. 

Recombination test 

To map the location and orientation of the M H V com¬ 
plementation groups, a modification of the approach 
previously described to map the location of the poliovi¬ 


rus and apthovirus ts mutants was developed (Cooper 
eta!., 1975; Lake etal., 1975; Mackenzie ef a/., 1975). 
Cultures of DBT cells were seeded at a density of 1 
x 1 0 6 cells/well (Falcon 6 well plate) in DMEM contain¬ 
ing 10% Nu-serum and 1% antibiotic/antimycotic. 
Combinations of ts mutants were mixed and inoculated 
onto DBT cells at a m.o.i. of 10 for each mutant. Virus 
was adsorbed for 1 hr at room temperature and the 
inoculum removed. Individual wells were washed 
gently two times with 2 ml of warm PBS to remove un¬ 
bound virus, and incubated at 32° for 16 hr in 2 ml of 
DMEM containing 10% Nu-serum and 1% antibiotic/ 
antimycotic. Virus progeny were harvested and frozen 
at -70°. Each cross was titered at both 32 and 39.5° by 
plaque assay and the recombination frequencies were 
calculated as the percentage of ts + virus present in the 
progeny from the following formula (Lake et ai, 1975; 
Mackenzie et ai, 1975): 

RF = (d g W r + 3)39.5 x i 0Q 

m 32 

(Afi) 3 9 . 5 was the titer of the cross at the nonpermissive 
temperature while (Afi) 32 represents the titer of the 
same cross at permissive temperature. (A + fi) 39 . 5 was 
the sum of the revertants of each parent strain assayed 
at the restrictive temperature. This particular approach 
allows recombination to occur without selective pres¬ 
sure at 32° and provides a more realistic assessment 
of the true recombination frequency by calculating the 
percentage of ts + virus present within a given popula¬ 
tion. The formula only measures single crossover 
events resulting in the ts + phenotype and does not ac¬ 
count for recombination events in the opposite direc¬ 
tion resulting in the double ts mutant phenotype. 

Recombination frequencies were compared to a ref¬ 
erence cross (LA7 x LA9) that was included in each 
experiment to obviate day to day variations. These mu¬ 
tants were chosen for study because they have low re¬ 
version frequencies and the mutations in LA7 and LA9 
have been mapped by T1 fingerprint analysis (Lai etal., 
1985; Makino et at., 1986b, 1987). Day to day varia¬ 
tions in recombination frequencies were standardized 
with the following formula: 


(Astd)(x/?exp) 




A std was the mean recombination frequency of the ref¬ 
erence cross LA7 x LA9 that was calculated from 12 
replicative crosses. xfl exp was the mean recombination 
frequency of the two ts mutants being standardized to 
the reference cross and xA eX p was the recombination 
frequency of the standard cross (LA7 X LA9) that had 
been calculated during this particular experiment. B st6 
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was the standardized recombination frequency of 
cross B. 

Construction of a MHV-A59 genetic 
recombination map 

The average recombination frequency for each cross 
was calculated by standardizing 5 to 10 individual 
crosses to the control cross (LA7 X LA9). Standard de¬ 
viations were calculated for each cross and the map 
locations positioned with respect to LA7 and LA9. Fol¬ 
lowing analysis of recombinant virus by ts + phenotype, 
the mutants were arranged according to standard ge¬ 
netic practice assuming (1) recombination events were 
proportional to the distance between alleles; (2) double 
crossover events were a rarity as compared to single 
crossover events. Previous data with apthovirus and 
poliovirus ts mutants support the validity of these as¬ 
sumptions (Cooper et al., 1975; Cooper, 1977; Lakeef 
a/., 1975; Mackenzie ef al., 1975). 

RESULTS 

Effect of the input multiplicity on the MHV-A59 
recombination frequency 

Previous studies with apthoviruses and polioviruses 
clearly demonstrated that the input multiplicity of infec¬ 
tion dramatically affects the final recombination fre¬ 
quency of a particular cross (Lakeef al., 1975; Macken¬ 
zie et al., 1975). These studies demonstrated that as 
the m.o.i. decreased, recombination frequencies de¬ 
creased accordingly. Polioviruses and apthoviruses 
cause cytopathic rounding of infected cells while MHV- 
A59 causes fusion and giant cell formation (Spaan et 
al., 1988). This feature of MHV-A59 infection could dra¬ 
matically alter the effect of input multiplicity on the final 
recombination frequency of a particular cross by allow¬ 
ing recombination between singly infected cells to oc¬ 
cur following fusion. To assess the effect of cell fusion 
and optimize the input multiplicity for maximum recom¬ 
bination frequencies, cells were inoculated with ts LA7 
and LA9 at varying multiplicities of infection and main¬ 
tained at 32°. Progeny were harvested at 16 hr postin¬ 
fection and titered at both the permissive and nonper- 
missive temperatures. In general, recombination fre¬ 
quencies were highest at m.o.i.s greater than 1.25 in 
which over 50% of the cells were co-infected with both 
ts mutants (Table 1). However, the most consistent re¬ 
sults were obtained with m.o.i.s greater than 5. In con¬ 
trast to results obtained with poliovirus and apthovi¬ 
ruses, recombination frequencies for MHV-A59 re¬ 
mained relatively high at input multiplicities in which 
only 7% of the cells were co-infected (m.o.i.—0.31) 
(Table 1). While direct evidence is lacking, the most 


TABLE 1 


Effects of Multiplicity of Infection on Intertypic MHV RNA 

Recombination Frequency 


m.o.i. 


Theoretical 
% of cells 
mixed ly infected 9 

Recombination 

frequency * 5 

ts LA9 

LA7 

20.00 

20.00 

100.00 

4.82 ±0.60 

10.00 

10.00 

100.00 

3.38 ±0.10 

5.00 

5.00 

98.70 

4.07 ± 1.71 

2.50 

2.50 

84.30 

3.26 ±0.90 

1.25 

1.25 

50.90 

4.28 ± 1.77 

0.62 

0.62 

21.30 

2.61 ± 1.60 

0.31 

0.31 

7.10 

1.53 ±0.60 

0.16 

0.16 

2.20 

0.81 ± 0.10 

0.08 

0.08 

0.60 

0.49 ±0.04 

0.04 

0.04 

0.15 

0.07 ± 0.02 

0.02 

0.02 

0.04 

ND C 

0.01 

0.01 

0.01 

ND 

20.00 

10.00 

100.00 

3.21 ±0.50 

20.00 

5.00 

99.30 

3.83 ±0.40 

20.00 

2.50 

91.80 

1.44 ±0.50 

20.00 

1.25 

71.30 

0.61 ± 0.01 

20.00 

0.62 

46.20 

0.61 ±0.05 

20.00 

0.10 

9.50 

0.05 ±0.02 


9 Calculated from the formula (1 - e Ma ) (1 - e Mb )100, where Ma 
and Mb are input multiplicities. 

'’Calculated from the formula rf - ^ 5>3 ™ 5 ~ / + 5>395 X 100; 

IAB) 32 

where AB 39 5 is titer at nonpermissive temperature, A -f 839.5 is the 
sum of the reversion frequencies, and AB^z is titer of the cross at 
permissive temperature. 
c Not detected. 


likely interpretation from these data is that coronavirus 
RNA recombination is a late event in the viral growth 
cycle which can occur after fusion of singly infected 
cells. 

Time course for the release of MHV-A59 
recombinant virus 

Previous studies with poliovirus indicate that the ap¬ 
pearance of recombinant ts + virus during a single 
growth cycle peaks early in infection and remains con¬ 
stant until the death of the cell (Lake etal., 1975; Mac¬ 
kenzie etal., 1975). These data suggest that poliovirus 
recombination occurs early in the virus growth cycle. It 
is unclear whether recombinant viruses are released 
early or late in the MHV-A59 infectious cycle. To ad¬ 
dress this question, recombination frequencies were 
calculated at different times in the virus growth cycle 
from duplicate samples of the cross LAI 6 x LAI 8, LA6 
X LAI 6, and NC2 X LA7. These mutants were chosen 
since they were of the RNA - phenotype and represent 
members in the same (LA6 and LAI 6—group A) or 
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Fig. 1. Growth curves during mixed-MHV infection. Cultures of 
cells were infected at a m.o.i. of f 0 each with two ts mutants repre¬ 
senting the same or different RNA" complementation group of MHV- 
A59. After 1 hr for adsorption, the monolayers were washed, media 
replaced, and incubated at 32° for 16-18 hr. Virus progeny were 
harvested at different times postinfection and titered by plaque assay 
at permissive or nonpermissive temperatures. Virus titers: LA6 
X LAI 6—32°, A; 39° A. LAI 6 X LAI 8— 32°, □; 39° ■ NC2 X LA7— 
32°, O; 39° •. 


different (LAI8—group E; NC2—group B; LA7— 
group F) complementation groups. The double mutant 
LAI 6 (group A/B) was used because of the low rever¬ 
sion frequency of this mutant at restrictive tempera¬ 
ture. Cultures of cells were co-infected with two ts mu¬ 
tants at a m.o.i. of 10 each and maintained at 32°. Virus 
titers were analyzed at different times postinfection at 
the permissive or restrictive temperatures. At 32°, all 
crosses replicated to virtually identical titers and were 
characterized by similar growth kinetics. However, the 
percentage of ts + virus was significantly higher be¬ 
tween mutants of different complementation groups 
(LAI 6 X LAI 8, and NC2 X LA7) (Fig. 1). The peak re¬ 
combination frequency was 5.05% between comple¬ 
mentation groups A/B and E (LAI 6 X LAI 8), 5.95% be¬ 
tween groups B and F (NC2 X LA7), and 0.2% between 
members of the same complementation group A (LA6 
x LAI 6). In contrast to results reported for poliovirus, 
coronavirus recombination frequencies increased 
throughout infection and peaked between 16 and 18 hr 
postinfection (Fig. 2). These data further suggest that 


MHV-A59 recombination may be a late event in the vi¬ 
rus growth cycle. 

Recombination versus reversion frequencies 

To assess the difference between the recombination 
versus reversion frequency, we performed a series of 
crosses between ts mutants representing different 
complementation groups of MHV-A59. The ratio of re- 
combinants/revertants ranged from 63 to 20,000 times 
higher than the sum of the spontaneous reversion fre¬ 
quency of each individual mutant (Table 2). The extent 
of the difference reflected both the stability of the mu¬ 
tants used in the cross and whether the mutants were 
members of the same or different complementation 
groups. Between different complementation groups, 
recombination frequencies between individual mu¬ 
tants were as much as 10-fold higher than those pre¬ 
viously reported for either polioviruses or apthoviruses, 
and support previous findings suggesting that high fre- 



Fig. 2 . Time course for the release of RNA recombinant viruses. 
Cultures of DBT cells were infected at a m.o.i. of 10 each with two 
different ts mutants as described in the legend to Fig. 2. RNA recom¬ 
bination frequencies were calculated by subtraction of reversion fre¬ 
quencies from singly infected cells as described under Methods. 
LA6 X LAI 6, □; LAI 6 X LAI 8, ■; NC2 X LA7, O. 
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TABLE 2 


Recombination vs Reversion Frequency during Mixed-MHV Infection 


Mutant 

32° 

39° 

Reversion frequency 3 

RF 6 

Proportion 0 

LA3 

9.0 X 10 7 

1.0 X 10 2 

i.i x io~ 6 

__ 

_ 

LA6 

8.5 X 10 7 

4.1 X 10 4 

4.8 X 10~ 4 

— 

— 

LA7 

5.5 X 10 7 

1.0 X 10 2 

1.8 X 10“ 6 

— 

— 

LA9 

8.6 X 10 7 

7.8 X 10 3 

9.0 X 10~ 4 

— 

— 

LA10 

5.6 X 10 7 

2.7 X 10 4 

2.7 X 10“ 4 

— 

— 

LAI 8 

1.0X 10’ 

1.9 X 10 6 

5.3 X 10' 6 

— 

— 

7X3 

4.4 X 10 7 

4.0 X 10 6 

— 

9.1 

20,000 

7X6 

9.6 X 10 7 

8.5 X 10 6 

— 

8.9 

207.0 

7X9 

9.7 X 10 7 

4.0 X 10 6 

— 

4.1 

512.0 

7X10 

6.5 X 10 7 

1.7 X 10 6 

— 

2.6 

63.0 

7X18 

2.9 X 10 7 

2.4 X 10 6 

— 

0.8 

2,182.0 


3 Calculated as 39732°. 

6 Recombination frequency calculated as described under Methods. 

0 Proportion is the number of ts" recombinant viruses divided by the number of revertant viruses titered at 39° from singly infected cells. 


quency recombination occurs during mixed MHV in¬ 
fection (Keck ef ai, 1987,1988b; Makinoefa/., 1986b). 

Obstacles to mapping the location of the MHV-A59 
complementation groups 

To obtain a mean recombination frequency for a 
standard cross, we chose two mutants whose genetic 
defects had been mapped by recombination studies. 
Complementation group E mutant LA7 maps 7 to 8KB 
from the 3' end of the genome in the S envelope glyco¬ 
protein gene while group C mutant LA9 maps at the 5' 
end of the genome in the putative polymerase gene (Lai 
ef a/., 1985; Makino ef ai, 1986b, 1987). The mean 
recombination frequency of the standard cross was 
obtained from 1 2 replicated crosses of LA7 x LA9 with 
a mean recombination frequency of 4.45 + 0.85 (Table 
3). To minimize day to day variations, crosses per¬ 
formed on different days included the standard cross 
(Table 3). 

Construction of a genetic recombination map for 
MHV-A59 complementation groups 

In contrast to studies with poliovirus and apthovirus 
ts mutants, MHV-A59 ts mutants are amenable to 
complementation analysis (Koolen ef ai, 1983; Lei- 
bowitz ef ai, 1982; Schaad ef ai, 1990) permitting the 
recognition of distinct gene functions and facilitating 
genetic mapping. Unique ts mutants, isolated from 
different mutagenesis experiments, were selected 
from each complementation group and crossed with 
the reference mutants LA7 or LA9. In addition, mutants 
from each complementation group were also crossed 
with mutants representing all other complementation 
groups. The recombination frequencies and standard 


deviations for each cross are summarized in Table 4. 
The distances between complementation groups are 
within statistical limits and permit the construction of a 
genetic map (Fig. 3). In support of findings presented 
in the accompanying paper, ts mutants from each 

TABLE 3 

Effect of Standardization on MHV Recombination Frequencies 


Standard 

LA7XLA9 

recombination Before Following 

frequency 3 standardization 6 standardization 0 


4.277 

Experiment 1 

LA7X LA9 

3.77 

4.45 

4.102 

LA7 X LA10 

2.30 

2.70 

4.448 

LA7XLA18 

0.81 

0.95 

5.714 

3.750 

3.980 

4.226 

Experiment 2 

LA7 X LA9 

4.24 

4.45 

LA7 X LAI 0 

1.56 

1.67 

3.425 

LA7XLA18 

1.46 

1.53 

5.700 

3.499 

6.025 

4.170 

4.445 ±0.845 





3 Observed mean of 12 replicated crosses. 

6 Observed mean of 3 replicated crosses. 

c Standardized against LA7 x LA9 by the formula 

(X^exp) 

= B sxdt where A std is the standard recombination frequency of LA7 
X LA9, xA exp is the recombination frequency of LA7 X LA9 cross in 
the experiment; xB exp is the recombination frequency of the cross 
prior to standardization, and £ std is the adjust recombination fre¬ 
quency. 
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TABLE 4 



Genetic Recombination Frequencies between Different MFIV Complementation Groups 

a 


Complementation 

groups 

LA3 

LA6 

LAI 6 

NC2 

LA9 

LAI 0 

LAI 8 

LA7 

A 









LA3 

— 

0.2 ± .02 

0.7 ±0.1 

1.3 ±0.3 

1.6 ±0.4 

5.5 ± 1.5 

5.6 ± 0.4 

8.5 ±1.4 

LA6 

— 

— 

1.1 ±0.2 

ND 

1.9 ±0.4 

5.8 ± 1.5 

4.2 ±0.8 

8.6 ±1.3 

LAI 6 

— 

— 

— 

ND 

1.3 ± 0.6 

2.7 ±0.3 

5.1 ± 1.5 

6.2 ±1.7 

B 









NC2 

— 

-— 

—- 

— 

1.5 ± 0.5 

3.8 ± 1.2 

4.3 ± 1.0 

6.4 ±1.8 

LAI 6 

— 

— 

— 

ND 

1.3 ± 0.6 

2.7 ±0.3 

5.1 ±1.5 

6.2 ±1.7 

C 









LA9 

— 

-— 

— 

— 

— 

1.4 ±0.6 

3.0 ± 1.8 

4.5 ±0.8 

D 









LAI 0 

— 

-— 

— 

— 

— 

—_ 

0.9 ±0.4 

2.6 ±0.9 

E 









LA18 

— 

— 

— 

— 

— 

— 

— 

1.21 ±0.3 

F 









LA7 

— 

— 

—— 

— 

— 

— 

— 

— 


a Recombination frequencies were calculated as described in the text. 


complementation group were found to map in similar 
regions of the genetic map. These data suggest that 
distinct cistrons were encoded in unique regions of the 
MHV genome. From the 5' end of the genome, the or¬ 
der of the MHV complementation groups was A, B, C, 
D, E, and F. Two small regions of overlap were de¬ 
tected between the group A/B and D/E mutants. Mu¬ 
tants within complementation group A mapped over a 
large domain suggesting that this gene product is very 
large. 

Crosses between ts mutants in the same comple¬ 
mentation group had characteristically low recombina¬ 
tion frequencies. Not surprisingly, higher recombina¬ 
tion frequencies were observed between mutants from 
different complementation groups and were as high as 
8.6%, or 17.2% assuming that crossover events occur 
in both directions. Since LA7 maps approximately 23- 
24 kb from the 5' end of the genome, a 1 % recombina¬ 
tion frequency occurs over 1300-1400 nucleotides of 
RNA. Assuming a constant recombination rate through 
the polymerase region, and that l_A6 maps near the 5' 
end of the genome, we can utilize the boundaries 
shown in the recombination map (Fig. 3) to predict the 
tentative nucleotide domains for each MHV comple¬ 
mentation group (Fig. 4). These data predict that the 
group A mutants map over a 4- to 6-kb stretch of RNA 
while groups B through E map over smaller domains. 
While the exact nucleotide domains cannot be pre¬ 
dicted from these data, it clearly demonstrates that ts 
mutants from the same complementation group map 


in similar nucleotide domains and indicates that all of 
the RNA" groups used in this study map in the polymer¬ 
ase gene at the 5' end of the genome. 

DISCUSSION 

Homologous recombination has been reported dur¬ 
ing poliovirus, apthoviruses, and coronavirus infection 
(Cooper eta/., 1975; Cooper, 1977; King era/., 1982, 
1985; Lai etal., 1985). For the former viruses, tempera¬ 
ture-sensitive mutations in the viral genome have been 
arranged into a genetic map by recombination analysis 
(Cooper, 1968; Cooper etal., 1975; Lake etal., 1975). 
In this report, we demonstrate that complementation 
groups of MHV-A59 can also be arranged into an addi¬ 
tive, linear, genetic map on the basis of their recombi¬ 
nation frequencies. This study represents the first re¬ 
port establishing a genetic recombination map for co- 
ronaviruses, and provides insight into the approximate 
number, size, and location of the virus genes which 
function in RNA synthesis. Moreover, utilizing well- 
characterized ts mutants, the total recombination fre¬ 
quency for the MHV genome can be predicted. 

The MHV genome is divided into seven or eight cod¬ 
ing regions (Fig. 4). Expression of each coding region 
is mediated through the transcription of genome- 
length or subgenomic mRNA which encode one or 
more proteins (Siddell, 1983; Spaan era/., 1988). For 
MHV, eight proteins have been identified which are en¬ 
coded in subgenomic mRNAs 2 through 7 (Armstrong 
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Fig. 3. Genetic recombination map for MHV-A59 complementa¬ 
tion groups A through E. Cultures of DBT cells were infected with 
different combinations of MHV-A59 ts mutants at a m.o.i. of 10 each, 
and recombination frequencies determined as described under 
Methods. Percentage recombination frequency is shown across the 
top of the figure. All mutants were initially positioned with respect to 
LA7 in the map and the recombination frequencies are shown by the 
solid arrows. Mutants were then crossed with LA9 and representa¬ 
tive mutants from each MHV-A59 complementation group and the 
boxed areas represent the predicted domains of each complementa¬ 
tion group. Hatched domains represent regions of overlap between 
the group A/B and D/E mutants. 


era/., 1984; Schmidt et al., 1987; Shieh era/., 1989; 
Skinner er a/., 1985; Skinner and Siddell, 1983, 1985; 
Spaan er a/., 1988). These include the three structural 
proteins N, M, and S which are encoded on mRNAs 7, 
6, and 3, respectively, and a 43-kDa glycoprotein (HE) 
encoded in the mRNA 2-1 of some strains of MHV. 
Four nonstructural proteins are encoded in the sub- 
genomic mRNAs. Viral mRNA2 encodes a 30-kDa non¬ 
structural protein while mRNA 4 is translated into a pro¬ 
tein of 15,000 Da molecular weight (pi 5). Two non¬ 
structural proteins designated pi 3 and plO are 
encoded in two open reading frames on mRNA 5. The 
function of these nonstructural genes is unknown. 

The gene(s) encoded at the 5' end of the genome is 
thought to function in regulating RNA synthesis (Spaan 
et al., 1988). The number of genetic functions encoded 
in the 23-kb polymerase region is unclear, but this re¬ 


gion contains the capacity to encode 700-800 kDa of 
protein (Pachukera/., 1989). In vitro, the MHV genomic 
RNA is translated into a 250-kDa precursor polyprotein 
which is subsequently cleaved into p28 and p220 pro¬ 
teins (Denison and Perlman, 1986,1987). The p28 pro¬ 
tein is encoded within the first 1.1 kb from the 5' end 
of the genome (Soe et at., 1987). An autoproteolytic 
activity, which cleaves the p28 protein from the p250 
kilodalton precursor, is encoded 3.9-5.3 kb from the 5' 
end of the genome (Baker et at., 1989). Size analysis 
predicts that the group A mutants map in a cistron 
which is much larger (4-6 kb) than the domain encod¬ 
ing the p28 protein (1.1 kb). 

Several lines of evidence indicate that the increase 
in ts + recombinant progeny that were scored in assays 
at the nonpermissive temperature represent true re¬ 
combinants and not aggregates or polyploid particles 
containing genomes with complementing mutants. 
First, recombination frequencies are additive and bio¬ 
chemical mapping of several mutants support the loca¬ 
tion of these cistrons (Keck et at., 1987, 1988c; Lai et 
al., 1985; Makinoefa/., 1986b, 1987). Second, the en¬ 
hancement in the ts + recombinant virus yields occur 
when the parental viruses have defects in the same cis¬ 
tron and are not capable of complementation (Schaad 
et al., 1990). Third, plaque-purified ts + progeny repli¬ 
cate efficiently at permissive and restrictive tempera¬ 
ture. 

The existence of a genetic map implies that cross¬ 
overs can occur throughout the MHV genome. Cross¬ 
over sites in recombinant virus have been demon¬ 
strated throughout the MHV genome (Keck et al., 
1987, 1988b,c; Lai et al., 1985; Makinoefa/., 1986b, 
1987). Intratypic recombination frequencies between 
group A and F mutants approach 8.6 or 17.2% assum¬ 
ing that a similar recombination frequency occurs in 
the opposite direction and results in the double ts phe¬ 
notype (Table 4). While the double ts mutant recombi¬ 
nant phenotype has not been identified, it seems likely 
that such recombinants exist. Temperature sensitive 
mutants mapping internally in the genetic map have ad¬ 
ditive recombination frequencies with ts mutants map¬ 
ping in either the 5' or 3' direction (Fig. 3). Thus, such 
mutants probably participate either as an acceptor or 
donor in template switching to produce recombinant 
virus. Similar arguments have also been used to predict 
the total recombination frequencies for the poliovirus 
and apthovirus genomes (Cooper ef al., 1975; King et 
al., 1987). Since LA7 maps approximately 24 kb from 
the 5' end of the 32-kb MHV genome and recombina¬ 
tion mapping data predict that complementation group 
A is too large to map in the p28 protein encoded within 
the first 1.1 kb from the 5' end of the genome (Figs. 3 
and 4), a 17.2% recombination frequency occurs over 
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an ~23 kb nucleotide domain or 1 %/1300 nucleotides 
of RNA. If we assume that the recombination frequency 
over the 5'-most 23 kb is equivalent to the recombina¬ 
tion frequency for the 3'-most 8 kb, these data predict 
that the recombination frequency for the entire ~32- 
kb MHV genome probably approaches 25% or more, 
and therefore would be the highest RNA recombination 
frequency described for a linear plus polarity RNA virus. 
It seems likely that the predicted recombination fre¬ 
quency probably represents an underestimate since it 
does not include the frequency of double crossover 
events (Keck etal., 1988c), the heterogeneity and dele¬ 
tion in the S glycoprotein gene (Parker et at., 1989), the 
possibility of hot spots of recombination, or crossover 
events resulting in noninfectious genomes. 

Genetic recombination maps between poliovirus ts 
mutants predict a total recombination frequency of 
2.2% over the 7433 nucleotide genome (Cooper et a!., 
1975; Cooper, 1977; Lake etal., 1975). Assuming that 
recombination events occur in both directions, these 
data predict that a 1 % recombination frequency occurs 
per 1700 nucleotide pairs of the poliovirus double- 
stranded RNA form (Cooper, 1977) and 1 %/1300 base 
pairs for MHV. Recombination frequencies for bacte¬ 
riophage T4 and Escherichia coli are 1 %/200 and 1 %/ 
1750 base pairs, respectively (Hayes, 1968). Thus, 
MHV high frequency RNA recombination is probably 
mediated by the large size of the genomic RNA and a 
difference in the RNA polymerase molecules of picor- 
naviruses and coronaviruses. This assumption is not 
surprising since the MHV polymerase probably tran¬ 
scribes leader RNA, dislodges from the template RNA, 


and reinitiates leader-primed transcription of subgeno¬ 
mic mRNAs from full-length or subgenomic minus 
strands (Baric etal., 1983, 1985; Makino etal., 1986a; 
Sawicki and Sawicki, 1990; Sethna et a!., 1989). It is 
unclear why high frequency RNA recombination oc¬ 
curs during mixed MHV infection, but it may provide a 
mechanism to circumvent potentially high RNA poly¬ 
merase error rates associated with RNA virus replica¬ 
tion (Holland et at., 1982) and a natural source of ge¬ 
netic diversity among coronaviruses (Parker et ai., 
1989). 

Reassortment frequencies for reovirus and influenza 
virus vary depending on the virus strain and genome 
segment, and approaches 2-33% and 42%, respec¬ 
tively (Fields, 1981; Scholtissek et ai, 1978). Thus, 
MHV high frequency RNA recombination approaches 
the reassortment frequency of segmented RNA vi¬ 
ruses. The mechanism for MHV high frequency recom¬ 
bination is unclear but may occur as a consequence of 
discontinuous, nonprocessive transcription during viral 
replication. If copy choice is the mechanism for MHV 
recombination, these data suggest that template 
switching occurs at high frequency. We have pre¬ 
viously demonstrated the presence of small RNA inter¬ 
mediates bound to or dissociated from the MHV repli¬ 
cative intermediate RNA (Baric et ai, 1983, 1985, 
1987). High frequency RNA recombination provides 
additional evidence that these small RNAs may repre¬ 
sent the functional intermediates of RNA transcription 
and recombination. 

Utilizing genetic recombination mapping tech¬ 
niques, we have formulated the first additive, linear, ge- 
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netic map of the MHV-A59 ts mutants. The orientation 
of these complementation groups from the 5' end of 
the genome was A, B. C, D, E, and F (Figs. 3 and 4). 
These data are based on two assumptions: (1) that the 
group A mutant LA6 maps near the 5' end of the ge¬ 
nome and (2) that the recombination frequency is rela¬ 
tively constant over the polymerase region. At this 
time, it is unclear whether the MHV recombination rate 
is linear over the length of the genome and this will 
affect the accuracy of the tentative map domains of the 
MHV complementation groups. However, the poliovi¬ 
rus and apthovirus recombination maps are linear and 
appear to represent the physical map locations of ts 
mutations (King era/., 1987). Complementation groups 
A, B, C, and D map in the 23-kb polymerase region at 
the 5' end of the genome. It also seems likely that the 
group E mutants map in the polymerase region since 
these mutants are incapable of transcribing detectable 
levels of either positive- or negative-stranded RNA at 
the restrictive temperature (Schaad et al., 1990) and 
map 5' to the mRNA 2 coding sequence (Luytjes et al., 
1988; Shieh et al., 1989). These data suggest that at 
least 13 genetic functions are encoded in the MHV ge¬ 
nome including five in the polymerase region, two each 
in mRNAs 2 and 5, and a single genetic function en¬ 
coded in mRNAs 3,4,6, and 7 (Fig. 4). Previous studies 
with sindbis virus ts mutants indicate that a single com¬ 
plementation group may encode a large multifunctional 
protein containing two or more functions (Hahn et al., 
1989; Mi eta!., 1989; Sawicki etai, 1981; Sawicki and 
Sawicki, 1985). In fact, the nsP2 protein of sindbis is 
probably involved in the initiation of 26 S RNA synthe¬ 
sis, acting as a protease that cleaves the nonstructural 
polyprotein precursors, and may be involved in the 
shutoff of minus-strand RNA synthesis (Sawicki and 
Sawicki, 1985; Hahn etai, 1989). Since the nucleotide 
domains of several MHV cistrons appear quite large 
(3-6 kb), and some mutations within a complementa¬ 
tion group map 1300-2600 base pairs apart (group 
A—LA3/LA6 and NC8/LA16), the MHV complementa¬ 
tion groups may also encode multiple functions. This 
possibility is currently under study. 

Previously we have shown that mutants in the group 
C cistron are probably blocked in the ability to tran¬ 
scribe negative strand RNA while the group D mutants 
are blocked in the ability to transcribe mRNA, but not 
leader RNA at restrictive temperature (Baric et al., 
1985; Schaad et at., 1990). Mutants from these com¬ 
plementation groups map over internal domains in the 
polymerase region. The recombination map also pre¬ 
dicts that the group A mutants map in the region of the 
autoproteolytic activity which is encoded ~3.9-5.3 kb 
from the 5' end of the genome (Baker et al., 1989). The 
function of the remaining complementation groups is 


presently under study. The location of a sixth RNA - 
complementation group identified by other investiga¬ 
tors is unclear, but could reside in the polymerase re¬ 
gion or internally in the viral genome (Koolen et al., 
1983; Leibowitz et al., 1982). By size analysis, the 
group A mutants map over a 4-6 kb nucleotide domain 
that could encode a protein of approximately 200 kDa 
molecular weight (Fig. 4). These data suggest that the 
group A mutants do not map within the p28 protein en¬ 
coded at the 5' end of the genome and predict that a 
sixth genetic function is encoded in this region. Alter¬ 
natively, an additional genetic function may map inter¬ 
nally in the polymerase gene since several large gaps 
of RNA are present within the genetic map (Fig. 4). 
Compton et al. (1987) have also demonstrated that 
anti-N monoclonal antibodies block in vitro transcrip¬ 
tion, and we have shown that N is tightly associated 
with the leader RNA and transcriptional complex (Baric 
etai., 1988; Stohlman ef a/., 1988). These data suggest 
that the N gene functions in RNA transcription and may 
represent the putative sixth RNA - complementation 
group not represented in our panel of ts mutants. 

The data presented in this report suggest that the 5' 
~23kb polymerase region of the coronavirus genome 
contains at least five separate polymerase gene func¬ 
tions that can be ablated by mutagenesis. Although the 
exact localization of these genetic functions may not 
be precise, the tentative map domains provide a mini¬ 
mal framework within which to begin characterization 
of the size, location, and function of the polymerase 
genes of MHV. 
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